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Abstract 

A geometric model for nonholonomic Lagrangian field theory is studied. The mul- 
tisymplectic approach to such a theory as well as the corresponding Cauchy formalism 
are discussed. It is shown that in both formulations, the relevant equations for the 
constrained system can be recovered by a suitable projection of the equations for the 
underlying free (i.e. unconstrained) Lagrangian system. 

1 Introduction 

During the past decades, much effort has been devoted to the differential geometric treat- 
ment of mechanical systems subjected to nonholonomic (i.e. velocity-dependent) constraints. 
To a large extent the growing interest in this field has been stimulated by its close connection 
to problems in control theory (see for instance the recent books [U IE] , which also contain 
extended bibliographies on the subject). As far as the study of Lagrangian (or Hamilto- 
nian) systems with nonholonomic constraints is concerned, one can essentially distinguish 
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between two different approaches. The first one, commonly called "nonholonomic mechan- 
ics", is based on d'Alembert's principle and on an additional rule that specifies a bundle 
of covectors along the constraint submanifold, representing the admissible reaction forces 
(or, alternatively, a subbundle of the tangent bundle, representing the admissible infinites- 
imal virtual displacements). The second one is a constrained variational approach, often 
coined "vakonomic mechanics". As is well-known, the dynamical equations generated by 
both approaches are in general not equivalent (for that matter, see for instance PQEHIS]). 

In this paper we will study an extension of nonholonomic mechanics to the treatment of 
(classical) Lagrangian field theories with external constraints. For the time being, we will 
confine ourselves to first-order field theories with constraints that depend on the independent 
variables (say the space-time coordinates), the fields and their first-order partial derivatives 
only. For the sake of clarity it should be emphasized that we use the term "classical field 
theory" here in a broad sense, i.e. it essentially refers to any physical (or other) system 
that can be described by the Euler-Lagrange equations derived from a Lagrangian density. 
The constraints involved are called 'external' to distinguish them from the type of con- 
straints appearing in gauge theories that find their origin in the possible degeneracy of the 
Lagrangian. 

Constrained field theories have already been studied extensively in the literature but, as far 
as the case of external constraints is concerned, those treatments are usually based on a field- 
theoretic analogue of vakonomic mechanics (see e.g. [21] )• The nonholonomic approach we 
are going to discuss here is a continuation of some work by E. Binz et al. [2 . We will refer to 
this theory as "nonholonomic field theory" mainly for its formal analogy with nonholonomic 
mechanics and also to distinguish it from the constrained variational approach. However, 
the term nonholonomic should perhaps be used here with some caution. Indeed, as can be 
inferred from a remark in [21], there does not seem to exist a general agreement yet on the 
precise meaning of the notion of nonholonomic constraint in field theory. In that respect, a 
recent paper by O. Krupkova ^1] may come to the rescue with a view on straightening out 
this matter, at least from a purely mathematical point of view. 

The mathematical framework for a nonholonomic field theory that has been proposed in 
[2j involves, among others, a generalization of d'Alembert's principle and of the so-called 
Chetaev rule that is commonly used in nonholonomic mechanics to characterize the bundle of 
constraint forms representing the admissible reaction forces. The constrained field equations, 
as well as an extension of the so-called De Donder-Weyl equations for classical field theories, 
are then derived in a finite-dimensional multisymplectic setting. The treatment in [2] also 
briefly deals with the Cauchy formalism for nonholonomic field theory. 

The purpose of the present paper is to continue and extend the work described in [2] . First of 
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all, we will relax the Chetaev-type conditions by allowing for constraint forms that need not 
be determined by the constraints themselves. Next, in the multisymplectic setting we will 
construct a kind of projection operator that maps solutions of the free (i.e. unconstrained) 
De Donder-Weyl equations into solutions of the modified (constrained) De Donder-Weyl 
equations that have been proposed in |2j. This is similar to the situation encountered in 
nonholonomic mechanics (see e.g. JH|). Finally, we also treat the Cauchy formalism for 
nonholonomic field theory. In particular, we will show that in case a global 'space-time' 
splitting of the base manifold can be fixed, the resulting structures and equations on the 
infinite-dimensional space of Cauchy data are reminiscent of those appearing in the geometric 
treatment of time-dependent nonholonomic mechanics. This is in full agreement with the 
results described for the unconstrained case by A. Santamana (see also |17j). 

The scheme of the paper is as follows. In the next section we recall some basic elements from 
the multisymplectic approach to (unconstrained) Lagrangian field theories on jet bundles, 
mainly in order to fix some of the notation that will be used. In Sections 3 and 4 we then 
discuss the construction of a nonholonomic model for first-order Lagrangian field theories 
with external constraints. The corresponding constrained field equations are given, as well 
as the modified De Donder-Weyl equations. Next, a projector is constructed that maps 
solutions of the De Donder-Weyl equations for the unconstrained Lagrangian system into 
solutions of the constrained De Donder-Weyl equations. In Section 5, as an example to 
illustrate the theory, we briefly consider the case of incompressible hydrodynamics. Section 
6 is devoted to the Cauchy formalism for nonholonomic field theory and in Section 7 we 
conclude with some general comments. 



2 Multisymplectic approach to Lagrangian field theo- 
ries 

There is an extended literature on the multisymplectic approach to classical field theories: 
see for instance El IT§] . where the interested reader can also find further references. 
Before reviewing some aspects of that theory, we first briefly recall some basic notions from 
jet bundle theory for which we mainly rely on the book by D.J. Saunders |27j . 

2.1 Jet bundles 

Consider a fibre bundle ir : Y — > X of rank m, whose base space X is assumed to be an 
oriented manifold of dimension n+ 1, equipped with a fixed volume form rj. (All manifolds, 
maps, vector fields and differential forms are assumed to be of class C°°). The first-order 
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jet bundle </% is the set of equivalence classes j^cfi consisting of those (local) sections of 
7i around a point x G X having the same Taylor expansion up to order one at x. J l ir is a 
(n + 1 + m + (n + l)m)-dimensional manifold, fibred over X with projection ix\ : J lr n — > X. 
In addition, J 1 it has the structure of an affine bundle over Y of rank (n + l)m, modelled on 
the vector bundle Vtt (g> ti*(T*X) — > Y, with Vn the bundle of 7r- vertical tangent vectors 
to Y, and with projection denoted by tc 10 : J 1 it — > Y. In particular, we have tti = 7r o 7r 10 . 

As far as coordinates on Y are concerned, we will always consider local bundle coordinates 
adapted to tt, written as y a ), fi = 1, . . . , n + 1; a — 1, . . . , m, and where the coordinate 
system (rr M ) on X is taken such that 77 = d n+l x := dx 1 A • • • A dx n+1 . The induced bundle 
coordinates on J 1 -k will then be denoted by y a , y'p. 

Given a vector field £ on Y, its first jet prolongation to J x ix will be written as In 
coordinates, if 



Q x fl Qya 

then 



dx^ dy a \dx^ dx>* J dy^ 

In terms of the volume form 77 on X, one can construct a special tensor field S v on J 1 ^, 
called the vertical endomorphism (see j2Z| for a formal definition). It is a vector- valued 
(n + l)-form that has the following expression in the coordinate system described above: 

5, = (dy a - y a v dx v ) A d% ~ ° 



where d n x il = i_d_d n+1 x. Note in passing that we will not make a notational distinction 



O.r 



between the volume form 7] on X and its pull-back to Y or to J x t\ under the respective 
projections. 

A 1-form 9 6 A 1 (J 1 7r) is said to be a contact 1-form whenever (j l <j))*9 = for each section 
(f) of 71. Locally, the module of contact forms is spanned by the 1-forms 9 a = dy a — y^ dx^ . 

In this paper, we will frequently deal with connections (in the sense of Ehresmann) on 71*1. 
Such a connection induces in particular a direct sum decomposition TJ x ix = Htti © Viz\, 
where Viri is the bundle of 7i"i-vertical tangent vectors to J 1 it and Htti is a complementary 
distribution, called the horizontal distribution. The corresponding horizontal projection op- 
erator will always be denoted by h. We recall that a connection on J 1 ^ uniquely determines 
a section T of the bundle (vri)i j0 : J x i\\ — > J 1 ^, called a jet field, and that, conversely, each 
jet field characterizes a unique connection on -K\ (see j2ZI for details). Therefore, in the sequel 
we will use the denominations 'connection' (on a jet bundle) and 'jet field' interchangeably. 
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An integral section of a jet field T : J l n — > J 1 ^ is a local section o of tt% such that 
j x a = To a. A jet field, or the associated connection, is called semi-holonomic if its integral 
sections a (if they exist!) are first jet prolongations of sections of n, i.e. a = j l cf) for some 
: X — > Y. Whenever a semi-holonomic jet field admits integral sections, it is called 
holonomic (or integrable). We now recall the conditions, in coordinates, for a jet field to be 
( semi- ) holonomic . 

In coordinates, the horizontal projector h of a connection on tti can be written as 

{ d d d \ 

h = dx» <g> U— + r a — + T a u — , (2) 

for some functions r^(x K , y b , y b ) and T'J ll/ (x K , y b , y^), and the associated jet field T then reads: 

T : JV — > jVa, {x^y a X) » (x^y a ,y a ^,T%). 

A local section a of 7r, with a(x) = (x M , <r a (x), a^(x)), is an integral section of the jet field 
T if 

da a da a u 

— = r* and — ^ = . 

From this expression, it easily follows that the connection will be semi-holonomic if = y® 
and holonomic (or integrable) if, in addition, r* = T° . 

We wish to emphasize that the previous discussion about connections on jet bundles can 
be presented in fully intrinsic terms: we refer again to for details. In particular, the 
condition for a connection (or jet field), with horizontal projector h, to be semi-holonomic, 
can be expressed by demanding that 

i\J) = 0, for each contact 1-form 9. (3) 

We will make use of this characterization later on. 



2.2 First-order Lagrangian field theories 

Given a fibre bundle ir : Y — > X and its associated first-order jet bundle J 1 ^, as considered 
above, we now briefly recall some aspects from the multisymplectic formulation of first-order 
Lagrangian field theory on J 1 ^, where the fields are the (local) sections of ir. Consider a 
Lagrangian density Lrj with L a smooth function (the Lagrangian) defined on J 1 -r . We say 
that L is regular if its Hessian matrix is non-degenerate, i.e. 

\dy-dy b J r 

at each point of J 1 ^. 
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Using the vertical endomorphism S v we can construct the following (n + l)-form on J Xr n: 

6 L = S*dL + Lrj 

and we then define the (n + 2)-form Ql = —dQ^, called the Poincare-Cartan form. The 
following coordinate expression for Ql will often be convenient: 

n L = ~g^dy a A d n+1 x - d l^—j A (dy a - y a v dx v ) A (4) 

If L is regular, which we will always assume in the sequel, the Poincare-Cartan form is a 
multisymplectic form according to the following definition. 

Definition 2.1 (see [(J I16L II 7| ) A closed m-form Q on a manifold M is called multisym- 
plectic if the mapping v x G T X M \— > i Vx Q(x) G A m_1 (T ;1 *M) is injective for all x G M. 

Note that symplectic forms (m = 2) and volume forms (m = dimM) are particular examples 
of multisymplectic forms and, moreover, these are the only two cases where the mappings in 
definition 12. II are surjective as well as injective (assuming M is finite dimensional). 

We now give a brief outline of the derivation of the Euler-Lagrange equations for first-order 
field theories. For more detailed treatments we refer to jSlCHllZZI- The action functional 
associated to the given Lagrangian density Lr\ is defined as 



5(0) = / (j^TLr] (5) 
u 



where U is an open subset of X with compact closure, and is a section of tt defined over 
U . A section is an extremal of (j3J) if 



= o, 

t=o 



Jt yHJt 

for any local flow {ft} on Y, consisting of diffeomorphisms defined on a neighborhood of 
4>{U), satisfying 7r o ip t = tt, ip = id., and keeping the boundary of 4>{U) fixed. Extremal 
sections are characterized by the following property. 

Proposition 2.2 A (local) section of tt is an extremal for the action J^p if and only if 

(j 1 0)*(^)=O, (6) 

for all vector fields ( on J 1 7T. 

Proof: See for instance B prop. 7.1.2], jTU thm. 3.1] or [Zj. \qed\ 
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In coordinates, equation (0) yields the well-known Euler-Lagrange equations: 

dL d ( dL 
dy a dx^ \dy^ 

A more general problem consists in looking for sections r of tti such that r*(^f2^) = for all 
vector fields ( on J 1 it (i.e. r need not be the prolongation of a section of 7r). This leads to the 
so-called De Donder equations of Lagrangian field theory. In case of a regular Lagrangian, 
the De Donder equations are equivalent to the Euler-Lagrange equations In this paper, 
we will mostly be concerned with a kind of linearized version of these equations that is more 
easy to handle and is obtained by looking for a connection on tti whose integral sections will 
be extremals of (0). More precisely we have the following important proposition: 

Proposition 2.3 Given a holonomic connection with horizontal projector h, then the inte- 
gral sections of the associated jet field are extremals of if and only if 

i h Q L = nVt L . (7) 

Proof: See j23 thm. 5.5.5] and [HJ . \qed\ 

Moreover, a simple coordinate computation shows that, given a regular Lagrangian L on 
J 1 it, a connection on tci satisfying (|7J) will automatically be semi-holonomic. Equation ((7j) 
is also referred to as the De Donder- Weyl equation of Lagrangian field theory. 



J = (a — 1, ... ,m). 



3 Nonholonomic Lagrangian field theory 

We now bring constraints into the picture. Suppose we have a Lagrangian system on J 1 ^, 
with regular Lagrangian L. Let C ■=— > J 1 ^ be a submanifold of J 1 it of codimension k, 
representing some external constraints imposed on the system. Although one can certainly 
consider more general situations (see e.g. [2]), for the sake of clarity, we will confine ourselves 
to the case that C projects onto the whole of Y, i.e. 711,0 (C) — Y, and that the restriction 
(7T lj0 )|c : C — > Y of 7Ti o to C is a (not necessarily affine) fibre bundle. In particular the 
latter is a quite restrictive condition but, with proper caution, one can probably carry out 
the further analysis under some weaker assumption. Finally, one could still require that C 
should not be the first jet bundle of a subbundle of ir, but since it does not really affect 
the present treatment we will not insist on that (i.e. our discussion also covers the case of 
"holonomic" constraints that do not essentially depend on the derivatives of the fields). 

Since C is a submanifold of J 1 7r, one may always find a covering of C consisting of open 
subsets U of J 1 ^, with U D C ^ 0, such that on each U G U there exist k functionally 
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independent smooth functions (p a that locally determine C, i.e. 

C n U = {7 G /Sr : ^(7) = for 1 < a < k}. 

We remark here that the assumption that (71*1 o)ic is a fibre bundle implies, in particular, 

dip 

that the matrix (— — )( r y) has maximal rank k at each point 7 £ C D Z7. 

With the given data — a Lagrangian and a constraint submanifold C — the question now 
arises how to construct a suitable constrained field theory. As pointed out in the introduc- 
tion, one possibility is to follow a constrained variational approach as, for instance, in [21]. 
Inspired by the situation in classical mechanics one may also think of another approach, 
called "nonholonomic" , which involves some additional ingredients (see For mechanical 
systems with nonholonomic constraints it is now well-known that the "vakonomic" and the 
"nonholonomic" equations of motion are not equivalent in general. Equivalence is achieved, 
however, if the constraint functions can be written as total time derivatives of velocity- 
independent constraints, i.e. if we are basically dealing with holonomic constraints (see e.g. 
[TU] and references therein). 

3.1 The constraint distribution 

Making a small digression to nonholonomic mechanics, we recall that the construction of 
the equations of motion of a mechanical system with nonholonomic constraints is based on 
the so-called d'Alembert principle and involves, among others, the specification of a suitable 
bundle of admissible "reaction forces" (and a corresponding bundle of admissible virtual 
velocities), defined along the constraint submanifold. This choice relies on an additional rule 
or principle. In nonholonomic mechanics it is quite common to use the so-called Chetaev 
principle, whereby the bundle of reaction forces is constructed directly in terms of the given 
constraints. In principle, however, the specification of the appropriate bundle of reaction 
forces (or virtual displacements), compatible with the given constraints, is problem depen- 
dent and need not necessarily be based on Chetaev's rule. For a critical discussion of this 
matter we refer to |22]; see a ls° EE]- 

Returning to the case of first-order field theory with external constraints, we now introduce a 
special subbundle F of rank k of the bundle of exterior (n + 1) -forms on J 1 ^ defined along the 
constraint submanifold C, where we recall that k is the co dimension of C. This bundle, which 
we will simply refer to as the bundle of constraint forms, will play a role similar to that of the 
bundle of reaction forces in nonholonomic mechanics. We define this bundle by considering 
a submodule T of rank k of the module of (n + l)-forms $, defined on a neighborhood of 
C, and which are n-horizontal and 1-contact, i.e. $ vanishes when contracted with any two 
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^-vertical vector fields and = for any section of it. In particular, one can find 

an open cover U of C such that on each open set U G U, the module T is generated by k 
independent (n + l)-forms $ Q that locally read 



= {CaT a {dy a - y a M) A d% = {C a yj a A d n x 



for some smooth functions (C a )a on U. Independence of the forms $ Q clearly implies that 
the (k x (n+ l)m)-matrix whose elements are the (C a )%, has constant maximal rank k. The 
bundle of constraint forms is then defined by 

F = |J F 7 with F 7 = {$(7) I $ G J"} . 

At this point, the reason for selecting a constraint bundle of the type described above is 
primarily based on the analogy with nonholonomic mechanics. 

Remark 3.1 In the authors have constructed the bundle of constraint forms by con- 
sidering a natural extension of the Chetaev-principle that is commonly used in mechanics 
when dealing with nonlinear nonholonomic constraints. More precisely, they define the local 
generators $ a of the bundle of constraint forms by putting 

$ Q := S*(d(p a ) , 

where the <p a are the local constraint functions (see the beginning of Section 3). One easily 
verifies that these $ a are indeed of the form with {C a )a — §pr- ^ n the case we are 
considering, the independence of these $ a is guaranteed by our initial assumption that C 
should have a fibre bundle structure over Y . 

As we will now show, the constraint bundle F gives rise to a distribution D along C, called 
the constraint distribution. As above, consider an open cover U of C such that on each 
U EU, the module T is generated by k independent {n + l)-forms $ a of the form (jSJ). 

Proposition 3.2 For each a, there exists a unique vector field ( a G X(U) such that 

i c n L = (9) 

Proof: Take Ca to be a tt^o- vertical vector field on U, i.e. Ca = (Ca)^a^- Herewith, equation 
reduces to 

(U « = (cx (10) 



9 



which determines the (Cq)^ uniquely, as L is supposed to be regular. This already proves 
the existence of a solution of Uniqueness then follows from the fact that is multi- 
symplectic. 

\qed\ 

The vector fields ( a span a fc-dimensional distribution D\j on U. It is not difficult to check 
that for any two open sets U,V G U with nonempty intersection, and for each 7 G U fl V, 
Du{l) — Dv{l)- Indeed, assume T is generated on U by k independent forms and on V 
by k independent forms $ a , then there exists a nonsingular matrix of functions r£ on £7 fl V 
such that 

$ a = rf $/3 . 

If we denote the corresponding generators of by C« and those of Dy by C a , ft readily 
follows from the previous proposition that 

Ca|c/nv = r cX/3[[/ny ' 

which proves that D v = Dy on [/ fl V. Consequently, the local distributions described in 
the previous proposition induce a well-defined (global) distribution D along the constraint 
submanifold C, whose sections are vr^o-vertical vector fields. Moreover, using a similar argu- 
ment as above one easily verifies that this distribution does not depend on the initial choice 
we made for an open cover U of C. 

3.2 The nonholonomic field equations 

Summarizing the above, we are looking for a field theory built on the following data: (i) a 
Lagrangian density Li] with regular Lagrangian L G C°°{J 1 tx)\ (ii) a constraint submanifold 
C C J 1 ^ that can be locally represented by equations of the form ip a (x fl ,y a ,y^) = 0, for 
a = l,...,k and where the matrix (d(p a /dy^) has maximal rank k; (iii) a bundle F of 
constraint forms and an induced constraint distribution D, both defined along C, whereby 
F is locally generated by k independent (n + l)-forms (JBJ), and D is defined according to the 
construction described in Proposition 13.21 

To complete our model for nonholonomic field theory, we now have to specify what the field 
equations are. Proceeding along the same lines as in [2] we introduce the following definition, 
using a generalization of d'Alembert's principle. 

Definition 3.3 A (local) section o of 71 : Y — > X, defined on an open set U C X with 
compact closure, is a solution of the constrained problem under consideration if j l a(U) C C 
and 

[ (j 1 a)*C^ ) L V = 0, 
Ju 
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for all ^-vertical vector fields £ on Y that vanish on the boundary of a(U) and such that 



%i)<I> = (*) 
for all sections $ of the bundle F of constraint forms. 

Putting £ = C, a (x,y)d/dy a and taking into account the expression (0) for the prolonged 
vector field it is easily seen that the condition (*) translates into 



{c a )T = o 



where the (C a )% are the coefficients of the constraint forms introduced in (JHJ). One can then 
verify that if o~(x) = (x^, o~ a (x)) is a solution of the constrained problem, then the functions 
a a (x) satisfy the following system of partial differential equations 



dL d ( dL 



dy a dx^ \dy^ 



K{C*T a (a = l,...,m), (n; 



Va (x^a a (x)^(x)) = (a = l,...,k). (12) 

As usual, the (a priori) unknown functions A^ play the role of 'Lagrangian multipliers'. The 
equations (fTTjl are called the nonholonomic field equations for the constrained problem. Note 
that if the bundle F of constraint forms is defined according to a Chetaev-type prescription 
(see Remark 3.1), then we recover the nonholonomic field equations derived in [2]. 

Let 1(F) be the ideal of differential forms, defined along C, generated by the constraint 
forms: i.e any element of T(F) is of the form £\ A, A for some G T and arbitrary 
differential forms Aj. Again proceeding along the same lines as in [2] we can formulate the 
following modification of the De Donder-Weyl problem for nonholonomic Lagrangian field 
theory: find a connection on tti : J lr K — > X with horizontal projector h such that along the 
constraint submanifold C 

i h Vl L - nVt L G 1(F) and ImhcTC. (13) 

For simplicity we will refer to ()13j) as the nonholonomic De Donder-Weyl equation. In 
coordinates, if we represent h by (J2J) one can easily check that the relation on the left of (|13|) 
leads to the following set of equations for the connection coefficients of the connection we 
are looking for: 



dy a dx T dyt T dy h dy a T TV dy b T dy a u Kv Uu) dy a dy b v r{a,a 
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for a = 1, . . . , m and /i = 1, . . . , n + 1 and some Lagrangian multipliers A". This should still 
be supplemented by the requirement that for any 7 £ C and any u £ T^J 1 ^, h(v) £ T 7 C. 
This is equivalent to requiring that h(v)({p a ) = for all v £ T c J 1 ^, where <p a (a — 1, . . . , k) 
are the (local) constraint functions. If, locally, h is written in the form (j2J, then the previous 
condition translates into the following additional equations for the connection coefficients in 
points of C: 

/i=l,.-.,n + l; a = l,...,fc. 

One can prove that in case of a regular Lagrangian, integral sections of a connection satisfying 
(f*H?|) will be 1-jet prolongations of solutions of the nonholonomic field equations (see [2] for 
details). 



4 The nonholonomic projector 

The purpose of the present section is to show that for a nonholonomic first-order field theory 
in the sense described above, one can construct, under an appropriate additional condition, 
a projection operator which maps solutions of the De Donder-Weyl equation (0) for the free 
(i.e. unconstrained) Lagrangian problem into solutions of the nonholonomic De Donder-Weyl 
equation (JT3|) . 

Given a constrained problem as described in the previous section, with regular Lagrangian 
L, constraint manifold C C J l n and constraint distribution D, we now impose the following 
compatibility condition: for each 7 £ C 

J D( 7 )nT 7 C = {0}. (14) 

If C is locally defined by k equations ip a (x tJ, ,y a ,y 1 ^) = and if D is locally generated by the 
vector fields ( a (see subsection I3.1J1 . a straightforward computation shows that the compat- 
ibility condition is satisfied iff 

det (C«(^)(t))^0, 

at each point 7 £ C. Indeed, take v £ T 7 C D D(j). Then v = v a ( a ('y), for some coefficients 
v a . On the other hand, = v((pp) = v a ( a ((p/3)(-y). Hence, if the matrix (C«(v 9 /?)(7)) is 
invertible, we may conclude that v — and the compatibility condition holds. The proof of 
the converse is similar. 

We now have the following result. 

Proposition 4.1 If the compatibility condition P^| ) holds, then at each point 7 £ C we have 
the decomposition 

T 7 J 1 7r = T 7 C© J D( 7 ). 
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Proof: The proof immediately follows from (|14|) and a simple counting of dimensions: 
dim T 7 C = dim T-yjV — k and dim .0(7) = k. \qed\ 

The direct sum decomposition of TcJ l n determines two complementary projection operators 
V and Q: 

V : T c J x tt -> TC and Q = I-V : T c J X n -> -D , 

where J is the identity on TcJ 1 ^. We will call P t/ie nonholonomic projector associated to 
the given constrained problem. 

Given a connection on tti such that the associated horizontal projector h is a solution of 
the free De Donder-Weyl equation (J7J), we will prove that the operator V o h|^ c ji w satisfies 
the constrained De Donder-Weyl equation Note that this operator is only defined 

along C and, therefore, strictly speaking it is not the horizontal projector of a connection on 
7Ti. However, one can show that its restriction to TC induces a genuine connection on the 
restricted bundle {it\)c '■ C — > X, and so the constrained De Donder-Weyl equation still 
makes sense for this kind of map. 

Lemma 4.2 The map V o hij^ji^ : TqJ 1 ^ — ► TC (c Tq J 1 ^), v i — > V(h(v)) is a projector 
whose restriction h-p to TC induces a connection on (^i)\c '■ C — > X. 

Proof: First of all, we check that for each 7 6 C the map V 1 o h 7 is a projector. Indeed, 
taking into account that ImQ = D is 71*1 ^-vertical, it follows that for all v G T 7 J l n 

(h 7 o Vry){v) = h 7 (ti) — (h 7 o Qj)(v) = h 7 (f ). 

and therefore 

(P 7 o h 7 )^ = P 7 o n 7 . 

The restriction hp of V o hi^ji^ to TC obviously is still a projector. The key point we now 
have to prove is that Im(h-p) is a complementary bundle to V{tz\)\c in TC, i.e. 

Im(hp)©V r ( 7 r 1 )| C = TC. (15) 

For that purpose we start by observing that along C we have TC fl Viti = V(tti)\c- m view 
of Proposition 14. II one can then easily derive the following direct sum decomposition: 

V(n 1 ) lc ®D = Vn 1 (along C). (16) 

Next, by taking into account the fact that the constraint distribution D is vertical, and 
therefore that hp(T 7 C) — (V o h)(T y J 1 7r) for every 7 G C, it is a routine exercise to verify 
that 

dim(V o hXT^Tr) = dimh^jV) . (17) 
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We now prove the direct sum decomposition (|15|). Take any v G TC with v G Im(h-p) D 
V{jti)\c, then there exists a vector u> G TC such that t> = V(h(w)) = h(w) — Q(h(w)). Since 
v is 7Ti-vertical, we conclude that h(w) = and, hence, v — 0. This already implies that 
Im(hp) n V(7Ti)|c = 0. The equality (fTK|) now follows from a simple dimensional argument. 
Indeed, relying on Proposition 14. II as well as on Qlfij) and (j!7j) . we have at each point 7 G C: 

dim(h-p(T 7 C)) + dimV^,(7Ti)| C = dim(h(T 7 J 1 ^)) + dim V^i — dim .0(7) 

= dim^jV) -dim£)(7) 
= dim T 7 C . 

This concludes the proof that hp = V o h|rc is the horizontal projector of a connection on 
(7Ti)| C . \qed\ 

Although (717 )|e : C — ► X is not a first-order jet bundle, we will say that a connection on 
( n i)\c, with associated horizontal projector h, is semi-holonomic if for each contact 1-form 
9 on J 1 ti 

10*6 = 0, (18) 

where j : C J 1 ^ is the canonical injection. Suppose r : X — > C is an integral section of 
a connection on (7Ti)ic, in the sense that Tt(T x X) C h(T T ( x .)C) for all x G Dom r. Then, if 
the given connection is semi-holonomic one can verify that, locally, r can be written as the 
first jet prolongation of a (local) section of 7r. 

As mentioned at the end of subsection 12.21 the regularity of L together with the fact that h 
satisfies the free De Donder-Weyl equation, imply that h is a semi-holonomic connection on 
J 1 ti. Herewith one can prove the following result. 

Lemma 4.3 The connection on (vri)|c defined in Lemma \4-2\ with horizontal projector hp, 
is semi-holonomic. 

Proof: We will use the fact that h is semi-holonomic and therefore satisfies ©■ Let 
t; E X,/ 1 ! be a 7170-vertical vector, then for any contact 1-form 9 on J 1 it we have that 
i v 9{^f) = 0. Now, for each v G TqJ 1 ^ we have that (V o h — h)(v) = — Q(h(v)) G D and, 
hence, (Voh — h)(v) is 7r 1)0 - vertical. Therefore i-poh@{v) = = for any contact 1-form 

9 and any v G TqJ 1 ^. From this one can readily deduce that h-p satisfies (I18J) and so we 
may conclude that the induced connection on {tt\)\c is indeed semi-holonomic. \qed\ 

We now arrive at the main result of this section. From now on, for ease of notation, we will 
use the projector V oh without further indication of its domain. The latter should be clear 
from the context. 



14 



Theorem 4.4 Consider a constrained problem of the type described above, with regular La- 
grangian L, constraint submanifold C C J l n and bundle of constraint forms F , and assume 
the compatibility condition \1J$ holds. Let h be the horizontal projector of a connection on 
tti, satisfying the free De Bonder- Weyl equation J7|) and let V be the nonholonomic projector 
associated to the constrained problem. Then, the projector V o h determines a solution of 
the constrained De Donder-Weyl problem M'J\) and restricts to the horizontal projector of a 
semi-holonomic connection on {tci)\c '■ C — > X. 

Proof: Along C we can rewrite the free De Donder-Weyl equation as 

Therefore, in order to prove that V o h satisfies the constrained De Donder-Weyl equation, 
we only need to verify that the right-hand side is an element of 1(F). 

We can write the projector h as h = dx^^H^ with = (d/dx^) -\-T°^(d / dy a ) -\-T a (d/dyf) 
(cf. (0))- Along C we can then put Q(H^) = X^C a for some functions and with the 
vector fields ( a as defined in Proposition 13.21 Then, at each point 7 G C and for any 
Vi, . . . , v n+ 2 € T 7 J 1 ^ we obtain 

n+2 

{iQ<A& L ){vx, . . .,v n+2 ) = o h)(^),wi, . . .,Vi, . . .,v n+2 ) 

i=\ 

n+2 

= ^(-l) m A^(^)(z Ca fi L )( Wl , . ,v n+2 ) 
i=i 

This shows that, along C, 

i Qob n L = -\«dx*A$ a ei:(F), 

which completes the proof of the first part of the theorem. 

The proof that Voh induces a semi-holonomic connection on (iri)\c follows from the previous 
lemmas 14.21 and 14.31 \qed\ 

Note that even in case a connection on 71*1, with horizontal projector h satisfying the free De 
Donder-Weyl equation, is holonomic (or integrable), the 'projected' semi-holonomic connec- 
tion h-p = V o h on (vri)c need not admit integral sections in general. 

5 An example: incompressible hydrodynamics 

As an example of a field theory with an external constraint, we consider the case of an 
incompressible fluid flow. This problem has already been treated for instance in [24J, us- 
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ing the constrained variational approach. From the point of view of "nonholonomic field 
theory" this is perhaps an a- typical example since, as we shall see, the constrained field 
equations resulting from the nonholonomic approach are essentially the same as those de- 
rived in [24J. The reason for this probably stems from the fact that, as we will show, the 
incompressibility constraint can be written as a divergence. Recall that for a mechanical 
system with a nonholonomic constraint that arises from a total time derivative of a function 
on the configuration space, the nonholonomic and the vakonomic equations of motion are 
equivalent. 

5.1 The constrained problem 

We will consider a simplified model of incompressible fluid motion in the sense that we will 
not bother about the technicalities related to precise domain and boundary conditions (see 
ffif] for a more detailed description of the geometric model). We identify the base space 
X with Kxl 3 e M 4 , with coordinates = {t,x l ) representing time t = x° and the 
material variables x i {i — 1, 2, 3). The volume form on X is the standard Euclidean volume 
1] := d 4 x = dt A dx l A dx 2 A dx 3 . For the total space we take Y = X x M 3 with coordinates 
(t, x l , y a ). The projection 7r : Y — ► X then reads ir(t, x l , y a ) = (t, x l ). In order to preserve 
some consistency with j2l|, we will denote the corresponding bundle coordinates on J 1 ^ by 
(t,x i ,y a ,v^,vf). 

Remark 5.1 In continuum mechanics it is common to denote the coordinates on X and 
Y by (tjX 1 ) and respectively, and the coordinates on J l ix by (t, X 1 , x l , v l , F}). 

However, we will not follow that convention here. 

In addition, we will equip the fibres M 3 of X (over the time axis) and of Y (over X) with the 
standard Euclidean metric, although one could replace them by more general Riemannian 
manifolds (see [231 121) • 

A section <p(t,x l ) = (t,x l ,(j) a (t,x)) of 7r can be seen as a map taking a material point x of 
the fluid and mapping it at each time t onto its position (fi a (t,x), (a = 1,2,3) in space. 
Following [23| we write the Lagrangian density as 

L{ 1 )d 4 x = - \\v \\ 2 pd 4 x-W(v?)pd A x, (19) 

where the function p = p(x) represents the material density, and W is the stored energy 
function. Note that W depends only on the v?, i.e. the "spatial" jet bundle coordinates. 
Next, we introduce the function J : J lr K — > K given by 

l 7( 7 )=det«( 7 )). 



16 



(note that (vf) is a square matrix). For any section of tt, J o j 1 ^ measures the volume 
change of a small fluid element under the 'flow' represented by <ft. In particular, the in- 
compressibility requirement can be expressed by the condition J7"(j'V) = 1' i- e - we have the 
constraint 

^( 7 ):=^(7)-l = 0, (20) 
defining the constraint submanifold C 

For the bundle F of constraint forms we take the line-bundle along C, generated by the 
4-form 

d>: = S* n {d V ) 

= ^(dy a -v a JxnAd 3 x^ 
= J{v- 1 )i(dy a -vtda^)A^x i) 
(i.e. we adopt the generalized Chetaev principle: see Remark [3. lj) . 



5.2 The nonholonomic field equations 

Before proceeding towards the field equations, we make the additional assumption that we 
are dealing with a barotropic fluid which, in particular, implies that W depends on the vf 
through J, i.e. W = W(J). The nonholonomic field equations (fTTj) for a barotropic fluid 
with Lagrangian (|19j) . subject to the incompressibility constraint (J2(Jj) and with constraint 
form $, become 

pSatj/o - (pW'Jiv-yj = \Jiv- 1 ): (a = 1,2,3) , 

which should be considered together with the constraint equation (p(t, x\ y a , v%, vf) = J — 1 = 

0. This should be compared with equation (4.8) in In that paper, the the field equations 
for an incompressible barotropic fluid were derived by means of a constrained variational 
approach. Since there is only one constraint equation, this approach gives rise to only 
one Lagrangian multiplier P, which can be interpreted as pressure. If we put Aj = J^, 
it is seen that, for the present example, the nonholonomic equations and the constrained 
variational equations are essentially the same. When thinking of the comparison between 
nonholonomic and vakonomic mechanics, the reason for this is to be found in the fact that 
the incompressibility constraint is determined by a divergence. More precisely, we have the 
following property. 

Proposition 5.2 The constraint function if can be written (locally) as a total divergence, 

1. e. there exist functions i/j 11 such that (p = 
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Proof: One can easily verify that 

d f dip \ dtp 



dx^ \ dv?, J dv° 



0. 



i.e. (p is a "null-Lagrangian" , which is equivalent to <p being a divergence (see e.g. [23 thm. 
4.7]). More directly, if we consider the functions 

tp° = and # = -Jy a (v~ 1 )i - x i , 

with v~ x the inverse of the matrix (vf), which are well-defined on a neighborhood of C, a 
rather tedious but straightforward computation shows that (p = dip^/dx 11 . \qed\ 

A detailed study of the comparison between the constrained variational approach and the 
nonholonomic approach to constrained field theories will be the subject of forthcoming work. 



5.3 The nonholonomic projector 

To illustrate some further concepts defined in the preceding sections, we now turn to the 
explicit form of the nonholonomic projector V for the example of incompressible fluid (not 
necessarily barotropic). As there is only one constraint, the bundle of constraint forces D 
is spanned by a single vector field ( = (®d/dv®. The coefficients of this vector field can be 
derived from (fTU|l where, in the present case, C% = dtp/dv^: 





If, for brevity, we denote the Hessian matrix of W with respect to the vf by 7i, then ( is the 
vector field along C given by 

Let us consider the function / := C(v 9 )> or explicitly 

f = (H~ 1 )%j\v- 1 y a (v-y b . 

For each 7 e C, f(j) 7^ from which it follows that the compatibility condition (|T4"|) holds. 
The nonholonomic projector V is then found to be 

v = i-r i d V ®(. 
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6 Cauchy formalism for nonholonomic field theory 



We will now describe the transition from the multisymplectic (covariant) treatment of non- 
holonomic field theory, discussed in the previous sections, to the formulation of the prob- 
lem on the space of Cauchy data. The Cauchy formalism for field theories is an infinite- 
dimensional analogue of classical dynamics of systems with a finite number of degrees of 
freedom. Instead of looking for sections of a bundle Y over an (n + l)-dimensional space- 
time manifold X (as in the covariant approach), one starts by introducing a space X of 
embeddings of a fixed 'Cauchy surface' into X. This space replaces the absolute time from 
Newtonian mechanics and, under suitable conditions, the system can then be described in 
terms of a particular vector field on an infinite-dimensional manifold Z, called the space of 
Cauchy data, which is a bundle over X. 

In our discussion of the Cauchy formalism for nonholonomic field theory, attention will be 
focussed on the case where the base manifold X admits a global splitting in 'space' and 'time'. 
It will be shown that, under appropriate assumptions, this Cauchy formalism reveals a close 
resemblance to the cosymplectic formulation of time-dependent nonholonomic mechanics (see 
e.g. This is in agreement with the results described in |26] for unconstrained Lagrangian 
field theory. Our aim is mainly to present the general idea, without entering into all technical 
details related to the geometry and analysis on infinite dimensional manifolds. 

Finally, it should be noted that only for hyperbolic partial differential equations it makes 
sense to consider initial value problems. In the remainder of this section we will therefore 
tacitly assume that the field equations we are dealing with, are hyperbolic in some suitable 
sense. We refer to jH] for a detailed analysis of this matter. 

6.1 The space of Cauchy data 

We first recall some basic aspects of the Cauchy formalism for Lagrangian field theories. We 
thereby closely follow the treatments presented in El EH] , to which we also refer for more 
details and further references on the subject. 

6.1.1 Generalities 

As before, we start from a fibre bundle it : Y — > X whose base space X is an (n + 1)- 
dimensional orientable manifold. Let M be an n-dimensional compact oriented manifold with 
volume form t]m- The pair (M,T]m) is called a Cauchy surface. A space of (parametrized) 
Cauchy surfaces X is then defined as a smooth manifold of embeddings t : M X. 
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Remark 6.1 Usually, X and M are taken to be manifolds with boundary and the embeddings 
r belonging to X are then assumed to map the interior, resp. boundary, of M into the 
interior, resp. boundary, of X. However, for the purpose of the present paper we will leave 
all considerations related to boundary aspects aside. 

In the sequel we will always assume, without loss of generality, that M has volume one, i.e. 

/ Vm = 1- (21) 

J M 

Points of M will usually be denoted by u. 

Given a space of Cauchy surfaces X, the space of Cauchy data Z is defined as a (infinite 
dimensional) manifold of embeddings from M into having the property that for each 
embedding k : M ■=— > J 1 7r, there exists a section of 7r and an element r of X such that 
k = j 1 (f)or. Finally, we define the space of Dirichlet data Y as consisting of those embeddings 
5 : M •— > Y having the property that there exists an element k of Z such that 5 = 7r li0 ° k. 

It is obvious from the previous definitions that the respective projections 7i~i : J x ix — > Y 
and ii : Y — > X induce the following natural projections: 

Z^Y^X. 

We further put ix i = tt o 7f 10 : the projection of Z onto X. 

The spaces X, Y, and Z can be equipped with the Whitney topology and can be made 
into (infinite-dimensional) manifolds (see 0). Because of the compactness of M, it does not 
matter whether one chooses the strong or the weak Whitney topology. However, as men- 
tioned at the beginning of this section, we will not unduly concern ourselves with technical 
issues related to the smooth nature of these manifolds and of the mappings and other objects 
defined on them. For a detailed discussion of the topology and differentiable structure on a 
space of differentiable mappings between manifolds, we refer to [T5] . 

The tangent space of X has a convenient geometrical interpretation (see [HJH2]): let r G 
X, then T T X can be identified with the space of sections of the pull-back bundle t*TX. 
Equivalently, a vector V T G T T X can be identified with a vector field along r, i.e. V T : M — > 
TX, u i — > V T {u) G T T ( U )X. Since r is a bijection (onto its image), one can still identify V T 
with a vector field on X, defined along r(M). Similar interpretations exist for elements of 
TY and TZ and will be freely used in the sequel. For instance, given k G Z, we will use the 
notation W K to indicate both an element of the tangent space T K Z and the corresponding 
vector field on J l n along k(M). 
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6.1.2 Existence of a global splitting of X 



For the further discussion we assume that there exists a global splitting of the base manifold 
X, induced by a diffeomorphism $ : K x M -* I. In physics, R is usually associated to 
time and M to (physical) space. An embedding r : M c — ► X is then called admissible if 
there exists a (necessarily unique) i G M such that t(u) = ty(t,u). We henceforth restrict 
X to be a manifold of admissible embeddings of M in X. This has the effect of reducing 
X to a 1-dimensional space diffeomorphic to M, with coordinate function denoted by t. The 
spaces of Cauchy data Z and of Dirichlet data Y, are then restricted accordingly. There is 
a canonically defined vector field H on X, given by 



S(r)(u) = 



for all u <E M, 

s=t 



where i is such that r(-) = •). In particular we have that (H, (it) = 1. The global time 
slicing of X allows us to consider a volume form 77 on X which, with some abuse of notation, 
can be written as 

77 = dt A tjm ■ 

(see e.g. [26 ). In the sequel we will always assume that X is oriented in terms of this volume 
form. 

An important property now is that, under the above assumptions (i.e. existence of a global 
splitting of X and X consisting of admissible embeddings only) one can prove that the space 
of Cauchy data Z is diffeomorphic to the jet bundle J 1 ^ (see |2T?1 par. 5.2.]). To make this 
diffeomorphism more explicit we note that there is a one-to-one correspondence between 
sections </> of 71*1 and sections (p of 7Ti: 

4>{x) = (p(r)(u), where x = \I/(£, u) and r(-) = ty(t, ■). 

We will frequently switch back and forth between both interpretations without warning, but 
we will stick to the notation "0" for a section of 7i"i and "y?" for the corresponding section 

of 7Ti. 



6.2 The unconstrained Lagrangian formalism 

Starting from a Lagrangian density Li] on J 1 ^, with regular Lagrangian L, the multisym- 
plectic (n + 2)-form induces a 2-form on the space Z of Cauchy data as follows. Let 
K E Z, and W K , W' K G T K Z, then put 

Q L (K)(H^,Wg := / K*(i W Jw>n L ), 
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where on the right-hand side, W K and W' K are interpreted as vector fields on J l n, defined 
along k(M) (see the end of subsection 6.1.1). Likewise, the (n + l)-form rj (pull-back of the 
volume form on X) induces a one-form fj on Z according to the prescription 

V(k)(W k ) := / K*(iw K v) 

J M 

for all k G Z, W K G T K Z. One can prove that both fl^ and fj are closed forms and, in 
particular, it turns out that CIl = — where O^ is the one-form on Z induced by (cf. 
[27)] for more details). 

As for the jet bundle 7%, one can show that the space of Cauchy data Z can be equipped 
with a 'vertical endomorphism' Sfj (see |2El section 5.2.3]). In the case under consideration, 
with X being 1-dimensional, Sfj is a vector valued one-form that can be defined as follows. 
Take any k G Z, with k = j 1 ^ ° T for some r G X and section <fi of 7T. In view of the 
identification between Z and J 1 ^ (see subsection 6.1.2), we can still represent k by jfy. For 
arbitrary W K G T K Z, we then put 

Sa(W K ) = {T jW n lfi (W K ) - T T ip o T^iW^Y , (22) 

where the superscript V denotes the natural vertical lift operation from TY into Vn^o. 
With the terminology used for vector fields on a first-order jet bundle, we will say that a 
vector field r on Z is a second- order vector field (shortly, a SODE) if 

Sfj{T) = and i r fj = 1 . (23) 



Consider a connection (or jet field) T on 7Ti : Y — > X, with horizontal projector h. One 
can then construct a vector field r on Z as follows. For k6 2, with k = j l cf) o r, define the 
vector T(k) G T k Z by 

r(«)(«) = h(TjV(SW(«))), (24) 

i.e. r(«)(u) G T K ( U ) J 1 ^ is the horizontal lift of S(r)(u) G T r ( u )X under the given connection 
T. We then have the following interesting property. 

Proposition 6.2 If T is a semi-holonomic connection on iti, then the vector field T on Z, 
defined by is a second-order vector field. 

Proof: For the contraction of F with fj we find that 

(i r fj) («) = / K*(ir( K )V) = / r *fe(r)^) = 1 , 
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where the last equality follows from the normalization assumption (|21j) and for the second 
equality we have used the fact that (with previous conventions) ir(K)V — tt* (i~.(T)v) and 
tti o k = t. Herewith, we have already shown that V verifies the second condition of (|2*3*j) . 

Next, we investigate the first condition of ()23j) . Since the given connection T is semi- 
holononomic, it is easily checked in coordinates that h satisfies 

T 7 7r li0 (h(w 7 )) = T 7 (0 o 7Ti)(f 7 ), (25) 

where 7 = j\§ and v 1 G T 7 J l ir. We now compute S^{T(k,)). With W K = T(k), the first term 
on the right-hand side of (}22|) becomes 

T k tti, (T(k))(u) = r K ( u )7ri, (r(K)(«)) 

= ^ K (n)7Tl,o(h(-ft' K ( u ))), 

where k = j^ip and where, for notational convenience, we have abbreviated Tj 1 0(S(t)(m)) 
by K K ( u y Using property (|23j) . we further obtain 

T K K 1;0 (r(K))(u) = T K(u) (0o 7ri)(ir K(M) ) 
= T k{u) <P(E(t)(u)), 

so that 

T^i,o(I>))=7>(~(r)), 

from which it follows that S^I^k)) = 0, which completes the proof that T defines a second- 
order ODE. 

\qed\ 

We then arrive at the following important result in the Cauchy formalism for (unconstrained) 
Lagrangian field theory. 

Theorem 6.3 Ifh satisfies the De Donder-Weyl equation Q), then the vector field T on Z, 
defined by fli^| ), satisfies the equations 

ir&L = and irfj = 1 . 
Proof: See |2H chapter 5]. \qed\ 

6.3 Nonholonomic constraints 

We now return to the nonholonomic setting described in sections H3 and 0] and, in particu- 
lar, we assume that the compatibility condition (fTlj) holds. In order to adapt the Cauchy 
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formalism, discussed in the previous subsection, to the nonholonomic case, we first define a 
subset C of Z as follows: 

C :— (k G Z | Ihik c c| . 

This set can be equipped with a smooth manifold structure such that C becomes a (infinite- 
dimensional) submanifold of Z. The tangent space to C at a point k is given by T K C = 
{W K G T K Z | W K {u) G T K(U) C}. 

For each k G C, let 

D K := G T K Z | Im^ K C , 

where .D is the constraint distribution along C. Putting 

d={Jd k 

one may verify that D determines a smooth distribution on Z along C. 

Next, for k G C and for each section a of the bundle F of constraint forms along C, we define 
an element a K of T*Z by 

a K (W K )= [ K*(i WK a), for all W K G T K Z. 

JM 

The set of all such covectors a K determines a subspace F K of T*Z and 
is a codistribution on Z along C. 

Since we assume that the given constrained problem satisfies the compatibility condition, 
we can use the nonholonomic projector V and the complementary projector Q = I — V (cf. 
Section 0} to define two operators V, Q : T^Z — > T^Z as follows. For each k G C and 
W K G T K Z, put 

V K (W K ) =VoW K (e Tj), Q k (W k ) = QoW K (e T K Z) . 

Using the properties of V and Q, it is not hard to check that, for each k G C, V k and Q K 
define complementary projectors in T K Z, i.e. 

(A) 2 = A, (Qk) 2 = Qk and V K + Q K = I K , 

with I K the identity on T K Z. This implies that T K Z = lmV K © lmQ K . Again relying on the 
definitions of C,D,V and Q, and on the properties of the nonholonomic projector V, one 
can prove that 

lmP K = T K C and ImQ K = D K . 
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Summarizing, we may conclude that under the given conditions we have the following de- 
composition of TZ along C: 

T d Z = TC®D. 

Let h be the horizontal projector of a connection T on 7Ti and let V denote the vector field 
on Z defined by ()24j) . The composition V o then determines a vector field on C, shortly 
denoted by V(T), and it is not difficult to see that it is precisely the vector field associated 
to the induced connection on (tii)\c with horizontal projector hp = ? o h (see Section HJ). 
We now have the following interesting result. 

Lemma 6.4 There exists a section a of F , such that 

if(r)&L = ir&L + 5- (26) 

Proof: For kgC and W K G T K Z, one can deduce from the definition of Ql that 

(ip, T) Q L )(K)(W K ) = / K*(i f{v){K) i WK VL L ). 

J M 

For the integrand on the right-hand side we have that, with u G M, 

1 V{V){k)(u) 1 W k (u)^L = «r(/t)(w)«W re («)^L - % Q(T){k)(u) 1 W k {u)Ql, 

where Q(r) is the vector field associated to Qoh (note that <2(r) is defined along C). Since 
Q(T)(k)(u) is an element of the constraint distribution D, the contraction with Ql yields a 
form a re ( u ) G F K ( u y Integration over M then gives (|2T)|) . \qed\ 

We have now collected all ingredients needed to formulate the main result of this section. 
Consider a constrained Lagrangian field theory, with regular Lagrangian L, with constraints 
verifying the appropriate conditions and such that the base manifold X admits a global 
space-time splitting. 

Theorem 6.5 Let h be a solution of the unconstrained De Donder-Weyl equation |7[) and 
let T be the corresponding second-order vector field on Z. Then, the vector field V(T) on C 
satisfies the following relations: 

ifirpL G F and V{Y) G TC. (27) 

Proof: If h satisfies the De Donder-Weyl equation, then the associated vector field T is 
contained in the kernel of Ql (see proposition I6.3|) . Expression (j2*o^) then proves the first 
part of (J27j) . The second part follows from the definition of C. \qed\ 

In addition, we note that V(T) is still a vector field of second-order type, due to propositions 
0and|O 
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To conclude, we have shown that under the appropriate assumptions, the Cauchy formalism 
for nonholonomic field theory leads to a vector field of 'second-order type' on the (infinite 
dimensional) subspace C of the space of Cauchy data Z, which can be written as a projection 
of the second-order vector field on Z associated to the free (unconstrained) Lagrangian 
system. 

7 Some final comments 

In this paper we have studied various aspects of nonholonomic Lagrangian field theory. 
Among others, we have shown that both in the multisymplectic approach and in the Cauchy 
formalism, the equations for the constrained system can be obtained by a projection of the 
equations for the original unconstrained Lagrangian system. 

While finalizing this paper we have come across a recent work of Olga Krupkova ([H]) in 
which nonholonomic Lagrangian field theory is discussed within the framework of a general 
study of partial differential equations with differential constraints. This paper — which 
differs both in purpose and methodology from ours — presents, among others, an interest- 
ing analysis of the various types of constraints that one may encounter when dealing with 
constrained exterior differential systems on fibred manifolds. 

The subject of nonholonomic field theory is still in full development. As far as the present 
study is concerned, there still remains some work to be done concerning the Cauchy formal- 
ism, mainly regarding the technicalities related to the infinite-dimensional manifold structure 
of the spaces Z and C. Another interesting matter that will be treated in future work, con- 
cerns the comparison between the constrained variational approach and the nonholonomic 
approach, i.e. the field theoretic analogue of the comparison between vakonomic and non- 
holonomic mechanics (see e.g. [TO]). 

Finally, an important challenge for future work will be the identification of some physically 
relevant examples to which nonholonomic field theory can be applied and for which, unlike 
the example of incompressible hydrodynamics treated in Section 5, the nonholonomic and the 
constrained variational approach are not "equivalent" . It is to be expected that interesting 
examples should come, for instance, from problems in elasticity (such as the rolling without 
slipping of a deformable body over a surface). 
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